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Abstract. We present a first combined analysis of the morphological and dynamical properties for the intermediate-mass 
Galaxy Evolution Sequence (IMAGES) sample. It is a representative sample of 52 z ~0.6 galaxies with M stc u from 1.5 to 15 
xlO lo M0 that possesses 3D resolved kinematics and HST deep imaging in at least two broad band filters. We aim at evaluat- 
ing the evolution of rotating spirals robustly since z ~0.6, as well as at testing the different schemes for classifying galaxies 
• i-H , morphologically. We used all the information provided by multi-band images, color maps, and 2D light fitting to assign a mor- 

phological class to each object. We divided our sample into spiral disks, peculiar objects, compact objects, and mergers. Using 
our morphological classification scheme, 4/5 of the identified spirals are rotating disks, and more than 4/5 of identified peculiar 
galaxies show complex kinematics, while automatic classification methods such as concentration-asymmetry and GINI-M20 
severely overestimate the fraction of relaxed disk galaxies. Using this methodology, we find that the fraction of undisturbed 
rotating spirals has increased by a factor ~ 2 during the past 6 Gyrs, a much higher fraction than was found previously based on 
morphologies alone. These rotating spiral disks are forming stars very rapidly, even doubling their stellar masses over the past 
6 Gyrs, while most of their stars were formed a few Gyrs earlier, which reveals a large gas supply. Because they are the likely 
progenitors of local spirals, we can conjecture how their properties are evolving. Their disks show some evidence of inside-out 
growth, and the gas supply/accretion is not random since the disk needs to be stable in order to match the local disk properties. 
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1. Introduction 

In the local universe, galaxy morphologies can be organized 
along the Hubble sequence. Both the physical and kinematic 
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properties of galaxies vary systematically with the Hubble type 
(Roberts & Haynes 1994). In this scheme, rotating disk galax- 
ies constitute the majority of the galaxy population. They rep- 
resent 70% of the intermediate-mass galaxy population, which 
itself includes at least 2/3 of the present-day stellar mass (see 
Ham mer et al. 20031 and references therein). 
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In the distant universe, morphological investigations 
based on HST imaging have brought observational ev- 
idence that a large fraction of galaxies have peculiar 
morphologies that do not fit into the standard Hubble se- 
quence (e.g. lAbr aham et al. 19961 Brinchma nn et al. 1 998 ; 



van den Bergh et al. 2001 Zheng et al. 2005 ). For instance, 
Zheng et al. (2005) conclude that at least 30% of the 
intermediate-mass galaxies at z=0.4-l have peculiar mor- 
phologies, whereas they represent less than few percent in the 
local universe. Over this fraction, more than 2/3 are luminous 
compact galaxies (LCGs - Ra wat et al. 20 07 ), an enigmatic 
population that has almost completely vanished in the local 
universe, with up to a factor ~10 decrease in comoving number 
density flWerk et al. 2004| >. 

While different observations suggest that most of the 
massive elliptical galaxies were in place prior to z=l 
(Jim enez et al. 200"71 IBernardi et al. 200 6). the formation and 
evolution of spiral galaxies is still under debate. Different 
physical processes can modify the galaxy properties over 
cosmic time, and the main hypothesis processes are: 

(i) the secular evolution with slow and continuous mat- 
ter accretion through the inter-galactic medium (e.g., 
ISemelin & Combes 20051 IBirnboim et al. 20071 1. 

(ii) minor mergers and accretion of small satellites 
dSomerville et al. 20011 1. 

(iii) more violent evolution through hierarchical merging (e.g., 
Hammer et al. 2005) as a function of lookback time. 

To help establish the relative importance of each process, 
both high-resolution imaging and integral field spectroscopy 
are required. Morphological studies can offer substantial clues 
to the merging rate dConselice et al. 20031 IBelletal. 2 006 
ILotz et al. 20061) . but kinematics studies using 3D spectroscopy 
appear to be a prerequisite to sampling the whole velocity 
field of individual galaxies and to directly distinguishing 
between interacting and non-interacting galaxies. The internal 
kinematics of galaxies are of primary interest because they 
probe the evolutionary state of these objects. 

In this context, we are pursuing an ESO Large Program us- 
ing the integral field capability of GIRAFFE (LP: IMAGES) 
to gather a complete and representative sample of velocity 
fields and dispersion maps of intermediate-mass galaxies at 
intermediate redshift. Galaxies are selected in different fields 
by their absolute J band magnitude (My(AB) < -20.3 - see 
Ravikumar et al. 20071 for more details). A first part of this 
sample was observed during the FLAMES/GIRAFFE guaran- 
teed time and is described in Flores et al. (2006), Puech et al. 
(2006a), Puech et al. (2006b), and Puech et al. (2007). The 
other part was observed in the frame of the large program 
IMAGES and is described in Yang et al. (2007) (hereafter Paper 
I). Based on this unique sample, Flores et al. (2006) and Paper 
I have developed a robust classification scheme to divide the 
sample into three distinct classes based on dynamical charac- 
teristics. Paper I finds that the whole sample can be distributed 
between these three kinematical classes because 43±12% of 
the galaxies have complex kinematics, 25+12% have perturbed 
rotation, and only 32+12% are consistent with pure rotation. 
Interestingly, objects supported by pure rotation all lie in the 



local Tully-Fisher relation dPuech et al. 2008al - hereafter Paper 
III), whereas objects identified as having perturbed or complex 
kinematics possess dynamical properties (specific angular mo- 
mentum, V/cr) that are statistically different from local disks 
and whose dispersions in properties are described by simula- 
tions of mergers dPuech et al. 20071) . In a first approach, Puech 
et al. (2007) suggest that galaxies with pure rotation could be 
associated with secular evolution, perturbed rotation with mi- 
nor mergers (Pue ch et al. 200 7b). and galaxies with complex 
kinematics could be the result of more violent events (e.g. ma- 
jor mergers). 

In the first part of this paper, we try to understand to what 
extent the morphological appearance of a galaxy reflects the 
underlying dynamical state. By combining the high-resolution, 
multi-band HST imaging with the 3D spectroscopy, we explore 
whether a correlation between morphology and kinematics ex- 
ists, like in the local universe. After describing the sample in 
Sect. 2, we present our derivation of a morphological classifi- 
cation of the sample in Sect. 3 and compare this classification 
with the kinematical classification in Sect. 4. In Sect 5 we com- 
pare the morphological method developed in this paper with 
automatic morphological classifiers. This morpho-kinematic al 
analysis will lead us to focus on one specific group: the rotat- 
ing spiral disks. In Sect. 6, we discuss the properties of these 
objects that dynamically and morphologically correspond to 
the local class of rotating spirals. In this paper we adopt the 
"Concordance" cosmological parameters of Hq = 70 km s _1 
Mpc~', Q.m = 0.3 and £2a = 0.7. Magnitudes are given in the 
AB-system. 

2. Data 

The whole sample is composed of 63 galaxies that were ob- 
served with the integral field units of FLAMES/GIRAFFE on 
the ESO-VLT. This data set comes from Paper I, and it com- 
bines the Flores et al. (2006) sample with the IMAGES large 
program sample, using the same selection criteria: Mj(AB) <- 
20.3 and EW([OII])>15A. The Flores et al. (2006) sam- 
ple includes 28 galaxies. Originally, spectra from the CFRS 
(Hammer et al. 1997) and from FORS1 and FORS2 data 
for the HDFS dVanzella et al. 20021 |Rigopoulou et al. 2005| > 
were used to select galaxies with the [Oil] emission lines de- 
tected. This first run of observations demonstrates that the 
FLAMES/GIRAFFE instrument was able to detect all galaxies 
with EW([OII])>15Aand l AB <23.5 after 8 to 13hrs of integra- 
tion time (see Flores et al. 2006 for more details). In the prepa- 
ration of the IMAGES large program, Ravikumar et al. (2007) 
obtained spectra of 580 galaxies in the GOODS-South field 
from VIMOS observations. From this catalog, combined with 
VVDS (Vimos VLT Deep Survey, ILe Fevre et al. 20041 ) and 
FORS2 ( Vanz ella et al. 20 06 ) spectroscopic data, Ravikumar et 
al. (2007) were able to draw a representative sample of 640 
galaxies at z <1 with Iab < 23.5. They show that, within 
the redshift range of 0.4< z <0.75, Iab < 23.5 galaxies in- 
clude almost all intermediate-mass galaxies (e.g. at least 95% 
of Mj(AB) <-20.3 galaxies). The 35 galaxies of paper I come 
from this sample. 

The Flores et al. (2006) sample and the IMAGES sample are 
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merged by applying the Mj(AB) <-20.3 and EW([OII])>15A 
selection criteria. This merged sample of 63 galaxies is repre- 
sentative of the luminosity function in the z=0.4-0.75 range, 
with a confidence level >99.99% (see Fig. 1 in Paper I). 
Because our morphological classification use color information 
(see Sect. 13.21 1. we chose to restrict ourselves to a smaller sam- 
ple (52 objects) with multi-band HST data (see Sect. I2.21 i. In 
the following we refer to the 63 galaxies as the "original sam- 
ple" and to the sub-sample of 52 galaxies having multi-band 
data as our "working sample". The working sample follows the 
z ~0.6 luminosity function with very similar confidence level. 

2.1. Kinematic observations 

A first part of the original sample (28 objects) was observed 
during FLAMES/GIRAFFE guaranteed time (ESO runs 07 1 .B- 
0322(A), 072.A-0169(A), and 75.B-0109(A)). The other part 
(35 objects) was observed in the frame of the large pro- 
gram IMAGES (ESO run 174.B-0328 - PI: F. Hammer). The 
[OH]/Lt3726,3729 doublet was used to derive both velocity 
fields and velocity dispersion maps. These maps were used to 
divide the sample into three distinct classes based on their dy- 
namical characteristics: 

- Rotating disks (RD): when the velocity field shows a ro- 
tation pattern that follows the optical major axis and the 
dispersion map show a peak near the dynamical center. 

- Perturbed rotators (PR): when the velocity field shows a 
rotation pattern but the dispersion map shows a peak not 
located at the dynamical center, or does not show any peak. 

- Complex kinematics (CK): when velocity field and disper- 
sion map show complex distributions, very different from 
that expected for simple rotating disks. 

A complete description of the methods and analyses for reduc- 
ing the data and classifying the galaxies is given in Flores et al. 
(2006) and Paper I. The final kinematic classification for each 
galaxy can be found in Table|4](col. 10). 

2.2. Complementary imaging data 

To carry out our morphological classification, we make use of 
all the HST data publicly available. Galaxies were selected in 
4 different fields, namely: the Canada France Redshift Survey 
03h and 22h, the Hubble Deep Field South and the Chandra 
Deep Field South. In the following, we describe each field and 
the corresponding HST data. 

2.2.1. Chandra Deep Field South. 

The CDFS (GOODS project) was observed with ACS in 
F435W (3 orbits), F606W (2.5 orbits), F775W (2.5 orbits), and 
F850W (5 orbits) bands (0.05"/pix - IGiavalisco et al. 20041 l. 
We used the publicly available version vl.O of the reduced, 
calibrated, stacked, mosaicked and drizzled images (drizzled 
pixel scale = 0.03"). Thirty-five objects of the original sample 
are located in the CDFS. All these galaxies are included in our 
working sample. 



2.2.2. Hubble Deep Field South. 

The HDFS was undertaken in 1998 (IWilliams et al. 20001 1. It 
consists of a deep principal field observed during 150 or- 
bits and 9 contiguous flanking fields observed for 2 orbit 
each dLucas et al. 20 03). For the principal field, WFPC2 imag- 
ing with the F300W, F450W, F606W, and F814W broad- 
band filters was performed (drizzled pixel scale = 0.04" - 
ICase rtano et al. 20001 l. The flanking fields were only observed 
with WFPC2 - F814W band (drizzled pixel scale = 0.05"/pix). 
Nine galaxies of the original sample are embedded in the 
HDFS: 5 are in the principal field, 2 are in FF2, and 2 in FF5. 
The 5 galaxies located in the principal field are included in the 
working sample. 



2.2.3. Canada-France Redshift Survey. 

The CFRS consist in 4 different fields located at 03h, lOh, 
14h, and 22h. It was partially observed with the WFPC2 
camera (0.1"/pix). A first campaign was in the 03h and 14h 
fields with F450W and F814W filters (ISchade etal. 1995k 
These observations were completed with F814W filter orbits, 
including the lOh and 22h field (Brinchma nn et al. 19 98). 
Finally, F606W and completion of F814W were carried out in 
the 03h and 14h fields ( [Zheng etal. 2004) . 
The original sample is divided between the 03h and 22h 
fields: 13 objects in the 03h field, 6 objects in the 22h one. As 
mentioned above, the CFRS was only partially observed by 
HST and only low-resolution, ground-based images (taken at 
CFHT, 0.2"/pix, |Schade et al. 1996t are available for 4 objects 
in the 22h field and 1 in the 03h. We do not include these 
objects in the working sample, as low-resolution data could 
lead to morphological misidentifications. The 2 remaining 
objects in the 22h field only have F814W images and are 
also not included in the working sample. The 12 remaining 
objects located in the 03h have either F606W and/or F450W 
images in addition to the F814W images and are included in 
the working sample (typical exposure time ~6400s). Finally, it 
is noteworthy that one object (CFRS031032) in the 03h field 
has been observed with ACS (F814W and F555W filters - 
0.05"/pix) during a gravitational lense program (Prop. 9744 - 
PLC. Kochanek). 

TableHJsummarizes the data available for the original sam- 
ple and the working sample. Middle row lists all the objects 
included in the original sample. Last row lists the objects with 
HST multi-band data included in the working sample. Despite 
the fact that the working sample includes galaxies observed at 
different resolutions, spatial samplings, and depths (ACS com- 
pared to WFPC2 and observed at different wavelengths), we 
find that none of these effects restricts our ability to homo- 
geneously classify all the galaxies morphologically. Moreover, 
note that the majority of the galaxies have high-resolution data 
taken with ACS with its better-sampled data and depth. 
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#of obj. 


03h 


22h 


HDFS 


CDFS 


Total 


Original Sample 


13 


6 


9 


35 


63 


Working Sample 


12 





5 


35 


52 



Table 1. Number of objects per field. From the 63 galaxies that eling of the galaxy flux assuming predefined light distributions, 
have been observed with GIRAFFE, we restrict ourselves to the We followed the same procedure as described in detail in Rawat 
52 objects with multi-band HST data. et al. (2007). Briefly, each galaxy is modeled as a combination 

of a bulge and a disk: the intensity profile of the bulge is mod- 
eled with a Sersic law and the disk as an exponential function. 
All the structural parameters (including Sersic index) are al- 
lowed to be free during the fitting process, except for the sky, 
which was held fixed at the estimate made by SExtractor. All 
neighboring objects within 5.0 arcsec of the target are simul- 
taneously fitted with a single Sersic component to avoid any 
flux contamination. Finally, during the fitting process, the star 
closest to each galaxy is used as the PSF for convolution (see 
Rawat et al. (2007) for more details). 

The results of the fitting procedure are a list of structural param- 
eters, an image of the modeled galaxy, and a residual map con- 
structed by the difference between real and modeled light distri- 
bution. This last map is used to compute a^- 2 that measures how 
far the modeled image is compared to the original light distri- 
bution. This parameter is generally used to estimate the quality 
of the fit. However, we find that this parameter alone is not suf- 
ficient for deciding whether a fit is reliable or not. To illustrate 
this point, Fig. [TJ shows the fitting procedure for two galaxies 
with their z-band images, fitting light models, and residuals. 
For these two objects, the software gives a similar level of con- 
fidence {x 1 = 1.38 and^- 2 = 1.46), whereas in the first case 
(left column) the x 1 is due to spiral arms and to real irregular- 
ities in the second (right column). In fact, the main limitation 
of parametric methods is that they use an "a priori" assumption 
on the light distribution, making them unsuitable for treating 
patchy light distribution often found in spiral galaxies, such as 
spiral arms. A spiral galaxy with prominent arms unavoidably 
produces a^ 2 > 1 because the spiral pattern is taken as if it was 
noise during the fitting. However, we can guage the appropri- 
ateness of the fit by carefully examining the residual maps in 
all of the high resolution images. A careful examination of each 
residual map allows us to define those that are relatively sym- 
metric due to spiral arm patterns that show the same spiral-like 
symmetry in all images (e.g., right column) from asymmetric 
ones (e.g. left column) that show asymmetries in all the images. 

Following this recommendation, we define a quality factor 
for each object to represent our confidence in the fitting result: 

Q=l: Secure fit: if the visual inspection of the fitting result 
and residual image shows that the Galfit solution is correct. 

Q=2: Possibly secure fit: if the visual inspection of the fit- 
ting result and residual image shows that the Galfit solution 
could be correct. 

Q=3: Unreliable fit or fit failed: if the visual inspec- 
tion shows that the Galfit type is not consistent with the 
original light distribution or if the fitting procedure crashes. 

We find that, for 40% of the galaxies in our sample, the 2D 
fit failed or seems unreliable (Q = 3) because light distribution 
presents overly distorted features. All these galaxies are mor- 
phologically peculiar objects. For galaxies that can be modeled 
by bulge+disk, we considered a segregation of galaxies into 
Galfit software ( Peng et al. 2002[ ), which performs a 2D mod- different Hubble types following IRawat et al. 20071 early type 



3. Morphological classifications 

We chose to follow a morphological analysis in three steps: (i) 
a surface brightness profile analysis was carried out to quantify 
structural parameters, (ii) we constructed a set of color maps 
for the whole sample, (iii) a morphological label was assigned 
to each object based on visual inspection of the images and 
color maps, and a detailed analysis of the structural parameters 
and physical properties derived from the first two steps. 



3.1. Surface brightness analysis 

The surface brightness profile analysis was essentially per- 
formed to derive the half-light radius and the bulge fraction 
(B/T). 



3.1.1. Half-light radius 

The half-light radius can be used as a compactness parameter. 
We consider that, under a half-light radius cutoff, objects are 
too compact to be clearly classified. These objects are fairly 
common at high z, and they represent a specific morphological 
category (Ra wat et al. 20 07). We chose to follow the compact- 
ness criteria defined in Melbourne et al. (2005), which classify 
as compact all objects with a half-light radius lower than 3kpc. 
Making a detailed morphological analysis for galaxies smaller 
than this size threshold, we found that our morphological 
analysis became unreliable. 

Two different techniques are used to derive the half-light 
radius. The first one is based on the IRAF task "Ellipse" 
(see Ha mmer et al. 20071 for a complete description of the 
procedure) and the second one on our own IDL procedure. In 
both cases, we adopted the method detailed in Bershady et al. 
(2000), which defines the total aperture as twice the Petrosian 
radius. Both methods give consistent results within an error of 
one pixel, which represents ~0.2kpc (~0.7kpc) at z=0.6 for 
ACS drizzled images (respectively WFPC2 in CFRS). 



3.1.2. Bulge fraction: B/T 

The B/T parameter is broadly correlated with the traditional 
Hubble type as early type galaxies have high B/T, whereas 
late types are expected to have low B/T. Then, it provides a 
valuable first guess for the morphological type of an object 
(Ken t 198Sb . To derive the B/T for each galaxy, we used the 
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Fig. 1. F850LP image, Galfit model image, residual image. 
For these two objects (J033237.54-274848.9 on the left and 
J033214. 97-275005. 5 on the right), the software gives a sim- 
ilar confidence level in the fitting process (% 2 -1.4). For 
J033237.54-274848.9 (left), the X 2 is high due to spiral arms, 
whereas it is due to the two clumps located southeast of the 
center for J033214.97-275005.5 (right). 



are for 0.8 < B/T<1, SO for 0.4 < B/T < 0.8 and late type for 
0.0 < B/T < 0.4. 

Both B/T and half-light radius are evaluated in each available 
filter. However, to minimize any effects of k-correction, we 
only keep those derived from F850LP (for CDFS) or F814W 
(for HDFS/CFRS) as it corresponds to the rest frame V band 
for all the objects, due to the narrow redshift range of the 
sample. Final results for the rest-frame V band are presented 
in Tableg](Cols. 4, 5, 6, & 9). 



3.2. Color maps 

We made use of the available multi-band images to con- 
struct two sets of color maps. The first set was obtained 
by combining, pixel by pixel, the magnitude in two ob- 
served bands (see Abra ham etal. 1999 Menant eau et al. 2 001 
Menante au et al. 20 04 Zheng et al. 2005 for examples). In the 



classification process, they are principally used to distinguish 
regular spirals with red bulges and blue disks (see Sect. 13.3b . 
Moreover, this color information can be compared with evolu- 
tionary synthesis models and then used to investigate the distri- 
bution of the stellar populations within the galaxies (see Sect. 



16.31 ). Physical properties, including whether some regions are 
dusty or starforming, can be derived from these color maps. 
This method provides a very direct way to identify the phys- 
ical nature of objects and compare distant galaxies with local 
galaxies in the Hubble sequence. The B435-Z850 color maps are 
produced for the 35 galaxies observed with ACS in the CDFS. 
For the objects located in the HDFS/CFRS, we obtained ei- 
ther B450-I814 or V606-I814 color maps. The pixels assigned to 
the color maps were derived following the method described in 
Zheng et al. (2004). This method allows quantitative measure- 
ments of the reliability of each pixel in the color map and then 
better constraints of the low signal-to-noise regions. Figure [T2l 
(second row) shows the color map images for the 52 galaxies 
with multi-band data available (see caption for details). 
As a complement to these quantitative color maps, we con- 
structed another set of "false three-color images" (see Figs. [3] 
and |4] for examples). These maps were constructed by com- 
bining three flux-weighted band images (B -V-z or B-V-I as 
blue-green-red). These maps cannot be used to derive physical 
parameters, however, in the morphological classification pro- 
cess, they are useful for deriving the global appearance of the 
galaxies. 

3.3. Visual classification 

A classification only based on B/T is difficult, principally 
because of the x 2 problem illustrated in Fig. [T] Moreover, 
it would not include the information provided by the color 
maps. To account for all the available information, we chose 
to visually classify all the galaxies. However, to reduce the 
subjectivity of such a task, we constructed a decision tree 
presented in Fig. [2] This decision tree was built to make the 
morphological classification simpler to understand. All the 
information provided by structural parameters, multi-band 
images, and color maps are included in the tree. At each step, a 
simple and unique criterion is considered. Following the arms 
of the tree step-by-step provides a reproductive and quantita- 
tive tool for performing the morphological classification. 

Galaxies were separated into four morphological classes as 
follows: 

1. Sp: the spiral disks, i.e., all objects with regular structures 
(arms) and a highly symmetric disk. Moreover, these spiral 
disk galaxies must show a bulge redder than the disk. 

2. Pec: the peculiar galaxies, i.e., objects with asymmetric fea- 
tures in the image or in the color map. As the peculiar class 
includes objects with different characteristics, we chose 
to divide this class into three subcategories: tadpole like 
(Pec/T), which are objects showing a knot at one end plus 
an extended tail ; suspected mergers (Pec/M), which are pe- 
culiar objects for which the irregularities could be associ- 
ated with merger/interaction events; and irregulars (Pec/Irr) 
for objects similar to local irregulars. 

3. C: the compact galaxies, i.e., all objects barely resolved and 
too concentrated to be decomposed into a bulge and a disk. 
We chose to classify as compact all objects with a half-light 
radius lower than 3kpc (see Sect. 13.1. il l. 
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4. M: the obvious merging/interacting systems, i.e., objects 
showing tidal tails, multiple cores, two components, etc. 
They generally also present very perturbed color maps, or 
extreme colors. 

To distribute our sample into the four classes, we followed 
the decision tree shown in Fig. [2] We first identified all the com- 
pact galaxies with the half-light radius information. Then, we 
investigated all the objects with B/T<0.4 to classify whether 
or not each galaxy is a disk. We made use of the color maps 
to classify the galaxies with a red bulge surrounded by a bluer 
symmetric disk. The remaining objects with B/T<0.4, but not 
satisfying the spiral disks criterion, were distributed between 
the Pec/Tadpoles and Pec/Irregular categories. To illustrate this 
process, Fig.[3]shows a three-color image of nine galaxies with 
B/T < 0.4 selected in the CDFS. The three galaxies in the first 
column satisfy the spiral disk criteria. The three galaxies in the 
second column clearly show more irregular features than those 
of the galaxies presented in the first column. These kinds of 
objects do not pass the "Spiral features?" test in the third row 
of the decision tree, and are then classified as Peculiar/Irregular 
galaxies. Similarly, objects like those in the third column do not 
pass the "Spiral features?" test, and as they show a knot plus an 
extended tail, they are assigned to the Peculiar/Tadpole galaxy 
class. 
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Fig. 3. B-V-z color images of nine galaxies selected in the 
CDFS sample with B/T<0.4. The three galaxies in the first 
column are classified as spiral disks, those in the second col- 
umn are Peculiar/Irregulars, and those in the third column are 
Peculiar/Tadpoles. 

Following the same procedure, we also investigated the 
0.4<B/T<0.8 range to look for regular SOs. We find no galax- 
ies satisfying this criterion. Finally, all the remaining objects 
for which B/T is unreliable or fit failed were distributed be- 
tween the peculiar and merger class. At this step, we identified 
the major mergers as galaxies with well-separated components 
with a strong color difference. Fig. [4] shows, for each morpho- 
logical class, a three-color image of a representative galaxy and 
Table[4](Col. 7) summarize the morphological classification as- 



signed to each object. 

Of course this method is not fully objective, and the distinc- 
tion between two classes may be ambiguous for some objects. 
The third object of the first column in Fig. [3] is an example of 
such a case, and it illustrates the limitations of our method. 
Nevertheless, these ambiguous objects only represent a low 
percentage of the whole sample, and the uncertainties are re- 
duced because 3 authors (BN, FH, AR) have independently 
classified the sample using the same set of criteria as outlined 
in Figure |2] 
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Fig. 4. B-V-z color images of six galaxies representative of 
the six morphological classes used in this study. From top 
to bottom and left to right: spiral disks, Peculiar/Irregular, 
Peculiar/Tadpole, Peculiar/Merger, compact galaxies, and ma- 
jor mergers. 



4. Comparing the morphological and kinematical 
classifications 

We compared the morphological classification derived in the 
previous section with the dynamical state of our objects as clas- 
sified in Paper I. As described in Sect. 12.11 galaxies were clas- 
sified into three kinematical class: RD, PR and CK (see Paper I 
for the kinematical classification). This kinematic classification 
was derived with a completely independent method, so it pro- 
vides a unique test of whether the morphological classification 
is consistent with the kinematic state of the galaxies. In Fig. [5] 
we report the comparison between the morphological classifi- 
cation derived with our method and the kinematic one. 

We find good agreement between morphological and kine- 
matical classifications: 80% of rotating disks are classified as 
spiral disks and more than 70% of complex kinematics are 
peculiar or merging galaxies. Among the galaxies classified 
as spiral disks, only one has complex kinematics and two 
have perturbed rotation. The spiral, whose kinematics were 
classified as complex (J033213. 06-274204. 8), is a galaxy for 
which the velocity field is not aligned with the optical axis. 
The two galaxies that show perturbed kinematics include one 
galaxy (J033248. 28-275028.9) that shows an elongated cr peak 
possibly due to a large bar and one (J033226. 23-274222. 8) 
that shows a sigma peak shifted from the dynamical center 
but located on a minor merger event (Pu ech et al. 200751) . In 
this last case, the morphology alone would be insufficient for 
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Fig. 2. Decision tree used to morphologically classify the sample. We assign each object a morphological label following step- 
by-step the arms of the tree. For each morphological class, we also report the number of galaxies. 



detecting such situations. Peculiar galaxies are mainly dis- 
tributed between PR and CK. Interestingly, galaxies identified 
as possible mergers or tadpoles are mainly CK (8/1 1), whereas 
Pec/Irregulars are mainly PR (9/13). Not surprisingly, all the 
galaxies identified as major merger are CK. Although the sam- 
ple is small, this may mean that our morphological classifica- 
tion correlates well with kinematical properties. Finally, com- 
pact galaxies are mostly CK (5/8) and only one is RD (see also 
IPuech et al. 2006b1 ). This last object is morphologically am- 
biguous as it possibly shows spiral arms, and has a very com- 
plex color distribution. 

To test our morphological classification against the kinematical 
one, we made a% 2 test between the two independent categories 
of rotating disks and complex kinematics. The y 2 test is used 
to determine the significance of differences between two inde- 
pendent groups. The hypothesis is that the two groups (RD and 
CK) differ with respect to the relative frequency with which 
their galaxies fall into different morphological classes. We find 
that the morphological classification of the two populations 
of RD and CK are inconsistent with a confidence level over 
99.98%. We also made the same test by considering the pecu- 
liar class as a single class, and we found a confidence level of 
99.99%. This suggests that our methodology is comparatively 
robust in classifying galaxies. 

5. Comparison with automatic classification 

We have shown in the previous section that it was possible to 
recover a good correlation between morphology and kinemat- 
ics following our classification scheme. More particularly, our 
method at least allows accurate separation of rotation disks 
from complex kinematics. Nevertheless, our method is not 



fully automatic and requires visual inspection of each galaxy. 
This human interaction unavoidably introduces subjectivity. 
In this section, we explore whether simple morphological 
parameters could efficiently distinguish RD from CK in a 
fully automatic way. As a first indication, we can use the 
B/T parameter derived by Galfit. For instance, we can as- 
sume a classification where all the objects with B/T<0.4 are 
considered as spirals. We find 28 galaxies satisfying such a 
criteria. However, beyond this number, only 14 are RD and 7 
are CK. We then conclude that a classification only based on 
B/T is not efficient at clearly distinguishing the RD from CK 
and would unavoidably lead to overestimating the number of 
rotating spirals. As illustrated in Fig. [3] a visual check of each 
object is necessary to clearly disentangle regular objects from 
more perturbed ones and to isolate only spiral disks. In fact, it 
appears that a full automatization of the decision tree would 
not be easy to produce. 

We also explore the "non-parametric" methods to find whether 
a correlation between these parameters and the kinematical 
state of a galaxy could be derived. To do so, we compute four 
parameters which are often used in morphological analyses: 
Concentration, Asymmetry, Gini and M20. The Concentration 
(C) parameter (Abraham et al. 19941 1 roughly correlates with 
the bulge over disk ratio (B/D) and the Asymmetry (A) 
parameter (Abrah am et al. 19 96) divides the sample between 
irregulars and more symmetric objects. More recently two 
other parameters have been introduced: the Gini coefficient 
(lAbra ham et a l. 20031 1 which is a measure of the relative 
distribution of galaxy pixel flux values, and the M20 parameter 
(ILotz et al. 2004K which is the relative second-order moment 
of the brightest 20% of a galaxy's pixels. 
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Fig. 5. Comparison between the morphological classification 
and the kinematical classification. Morphological classes are 
quoted in the bottom of X-coordinate, kinematical classes in 
Y-coordinate. RD refers to rotating disks, PR to perturbed rota- 
tions and CK to complex kinematics. 



For C and A, pixels assigned to each galaxy are defined at 1 .5<x 
above background. C and A are measured following Abraham 
et al. (1996) method. Given the narrow redshift range in the 
sample, we do not attempt to correct the concentration as 
described in Brinchmann et al. (1998) (Appendix A). Gini and 
M20 are measured following Lotz et al. (2004). 

In Fig.[6]we show how the sample is distributed in a log(A)- 
log(C) plane (top) and in a Gini-M20 plane (bottom). We find 
that the segregation between different morphological classes is 
partially consistent with our classification. For instance, there 
is a clear offset between the compact galaxies (black diamonds) 
and the spiral disks (blue circles), whereas the distinction be- 
tween peculiar (green squares) and spiral disks is less obvious. 
All the spiral disks indeed lie in specific regions in both dia- 
grams. However, several objects that have been identified by 
our method as peculiar (or even merger (red triangles)) galax- 
ies are mixed with the spiral disks. These last peculiar galaxies 
indeed show symmetric features or smooth light distribution, 
which explain why they lie near spiral disks. However, as illus- 
trated in Fig. [3] they are clearly inconsistent with our spiral disk 
criterion. It could be then difficult to define a criteria to isolate 
the spiral disks from other morphological types. As an exam- 
ple, we add the limits defined in Abraham et al. (1996) and 
Lotz et al. (2008) in both diagrams. These limits have been cal- 
ibrated on local samples and are generally used to distinguish 
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Fig. 6. top: log(A)-log(C) diagram. Bottom: Gini-M20 dia- 
gram. Different symbols refer to the morphological classes: 
spiral disks are represented as blue circles, peculiars as green 
squares, compact galaxies as black diamonds, and mergers as 
red triangles. 



three or more rough morphological types. With these limits, 
we find that the agreement for disks is good, since most of the 
objects designed as spiral disks are lying in the spiral (respec- 
tively Sb/bc) domain of both diagrams. The agreement is also 
good for the compact galaxies, as they almost all lie in or near 
the early type domain, which is consistent with the findings of 
Abraham et al. (1996). Note, however, that these objects could 
not be directly classified as early type objects because they 
show a relatively complex and/or blue color distribution (see 
sect. 16. 3 1 and Fig.fT2l. For peculiar objects, the assignments are 
not very accurate, and around half of these objects they do not 
lie in the peculiar region: Pec/Irr/Merg for C- A and Sd/d/Irr and 
Merger for Gini-M20. These results are consistent with those 
of Cassata et al. (2005), who find that these methods (even 
if efficient in disentangling early-type from late-type) gener- 
ally failed to resolve the different classes contributing to late- 
type galaxies (see also Conse lice et al. 2 003 Yagi et al. 2006 
Ivan derWel2008l) . 

To test the impact of this limitation on the morpho- 



kinematical correlation, we assume a morphological classifi- 
cation based only on the boundaries defined in Abraham et al. 
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(1996) and Lotz et al. (2008). All the objects lying in the spiral 
region for the log(A)-log(C) plane (respectively Sb/bc for the 
Gini-M20 plane) are assumed to be spirals, and all objects lying 
in the Pec/Irr/Merg (respectively Sd/d/lrr and Merger for Gini- 
M20) are assumed to be peculiar. The compact galaxies are 
included independently as they are identified when using the 
half-light radius. All the remaining galaxies lying in the E/SO 
domain are not included, because the working sample does not 
contain early type galaxies. In Fig.|7]we report the comparison 
between these classifications and the kinematical one. 
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Fig. 7. Comparison between the concentration-assymetry (dot- 
ted filled histogram), Gini-M20 (solid histogram) mor- 
phological classification, and the kinematical classification. 
Morphological classes are quoted in the bottom of X- 
coordinate, kinematical classes in Y-coordinate. RD refers to 
rotating disks, PR to perturbed rotations and CK to complex 
kinematics. The line for the CA and Gini/M20 exactly overlap 
for the RD types. 

In Fig. |7l it appears that the morpho-kinematical correla- 
tion is less obvious with the automatic methods. The same x 1 
test as in Sect.|4]gives 50% and 90% confidence for the log(A)- 
log(C) and Gini-M20 planes, respectively. For comparison, our 
classification considering the peculiars and mergers as a sin- 
gle peculiar class gives 99.9%. Whereas these methods almost 
recover the number of rotating spirals correctly, they mix the 
spiral disks with some peculiar galaxies: several galaxies lying 
in the spiral region have complex kinematics. This result in- 
deed confirms the finding of Conselice et al. (2003), Conselice 
et al. (2005), and more recently Atkinson et al. (2007), and 
Kassin et al. (2007). For example, Conselice et al. (2003) find 
that the asymmetry index is not sensitive to all phases of the 
merging process, particularly those in the beginning or at the 



end of the mergers. They argue that the total number of merg- 
ers would be underestimated by a factor of two, which is sim- 
ilar to what we find with our analysis. Kassin et al. (2007) re- 
port that some galaxies that were fairly normal according to 
Gini/M20 looked more disturbed or compact by visual inspec- 
tion. These objects generally present larger scatter in the Tully- 
Fisher relation (TFR) than visually identified disks. Combined 
with the result of Flores et al. (2006) showing that perturbed 
and complex kinematics are driving the scatter of the TFR, this 
also suggests that automatic methods mix the spiral disks with 
complex kinematics. Similar to B/T, a method only based on 
a few morphological parameters will unavoidably overestimate 
the number of rotating spirals, and it fails to uniquely identify 
all the CKs. 

6. Rotating spiral disks 

Our classification scheme is efficient in isolating different 
morpho-kinematical classes. We therefore use this method to 
focus on a particular subsample defined by the rotating spiral 
disks. This subsample includes all the galaxies sharing the spi- 
ral disk morphology and a relatively relaxed dynamical state 
(rotating disks and perturbed rotations). In doing so we have 
assumed that a disk, even with a slight kinematic peculiar struc- 
ture (i.e. PR), may be considered as a disk. This includes two 
galaxies, J033226. 23-274222. 8 (a minor merger, see Puech et 
al, 2007) and J033248.28-275028.9. Over the 52 galaxies of 
our sample, we find 14 objects sharing the spiral disk morphol- 
ogy and a rotating velocity field. These galaxies are the 14 first 
objects of Table|4](see also Tableland Fig.fTZli. they are partic- 
ularly interesting because they can be directly compared with 
the local rotating spirals. 

6.1. Sizes and B/T 

In terms of size and B/T distributions, these galaxies are con- 
sistent with the local late type spirals. The disk scale lengths 
derived from Galfit range from 2.5kpc to 6.5kpc with a mean 
value of 4.1kpc. For comparison, the disk scale length of the 
Milky Way is estimated to be between 2kpc and 3kpc, whereas 
the disk scale length of M3 1 ranges between 5kpc and 6kpc 
(see Sect. 2.1 in Hammer et al. 20071 and references therein). 
This indicates that the rotating spiral disks are roughly in the 
same size range compared to the local spirals. In terms of B/T, 
all these rotating spirals are distributed between 0.01 and 0.17. 
Around half of these objects have B/T<0. 1 consistent with Sc 
type or later, and the other half are in the range 0. 1 <B/T<0.2, 
consistent with Sb type. 

6.2. Fractional number of rotating spiral disks 

Over the 52 galaxies of the working sample, we find 14 
objects (27%) with both a spiral disk morphology and a 
simple rotational velocity field. As discussed in Sect. [2j 
the sample is representative of intermediate-mass galax- 
ies with EW([OII])>15A at intermediate redshift. Let us 
now consider the general population of galaxies at z=0.4 
- 0.8, with M ste u > 1.5xlO lo M , including those with 
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EW([OII])<15A. Following Hammer et al. (1997), the frac- 
tion of intermediate-mass galaxies with EW([OII])>15A rep- 
resents 60% of intermediate-mass galaxies. Then, the fraction 
of rotating spiral disks depends on how we distribute the re- 
maining 40% of those galaxies with EW([OII])<15A between 
the different morphological classes. To do so, we first assume 
that these quiescent galaxies are mainly E/SO and evolved 
spirals (see for instance Table 1. in ICassata et al. 20 05 ). We 
then assume a population of 23% of E/S0 consistent with 
a similar morphological study based on a sample of 111, 
0.4< z <1 intermediate-mass galaxies (Zheng et al. 2005). 
This fraction is also consistent with other several stud- 
ies (e.g. van den Bergh et al. 200 1| Brinchma nn et al. 1 998 
ICassata et al. 20051 IConselice et al. 20051 ILotz et al. 20081 . 
This leads to a fraction of 17% of evolved spirals with 
EW([OII])<15A. Accounting for these corrections, we find that 
33% (= 27x0.6 + 17)% of intermediate-mass galaxies are ro- 
tating spiral disks at z ~0.6. In Table [2] we report the distribu- 
tion of all the morpho-kinematical classes. RspD is for rotating 
spiral disks, C for compact galaxies and M for mergers. The 
morphological distribution for galaxies with EW([OII])<15A 
is from (Zheng et al. 2005). 



Table 2. 

classes. 



Fraction number of the different morphological 
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23% 17% 
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6.3. Colors 

Based on the color information provided by the color maps, 
we now investigate the distribution of the stellar populations 
within each galaxy. To do so, we compare the observed col- 
ors of our objects with modeled color-redshift evolution curves 
as predicted by Bruzual-Charlot stellar population synthesis 
code GALAXEV dBruzual & Chariot, 2003] ) . We can compute 
3 models of galaxies providing B435-zggo (for galaxies located 
in the CDFS) and V 606 -l8i4 (for CFRS and HDFS) observed 
colors at different redshift. The 3 spectral template are (i) an 
instantaneous burst corresponding to elliptical type objects; (ii) 
an exponentially decaying star formation rate (SFR) with e- 
folding time scale of lGyr, corresponding to SO like objects; 
(iii) an exponentially decaying SFR with e-folding time scale 
of 7Gyr, corresponding to late type like object. All these mod- 
els assume a solar metallicity and an epoch of formation at z=5. 
We also compute the theoretical color of a young starburst as- 
suming a power law spectrum L v oc v _1 and add the evolution 
curve of an Sbc galaxy from Benitez et al. (2004) templates. 
Figure|8]shows these evolution curves, as well as the integrated 
color of each galaxy. For comparison, we include all the mor- 
phological classes. 

In Fig. [8] almost all our objects lie between the starburst 
and the Sbc evolution model. No specific trends between differ- 
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Fig. 8. B435-Z850 (top) and V606-Isi4 (bottom) observed color 
as a function of redshift. Three solar metallicity, beginning 
their evolution at z-5 galaxy models, including an elliptical 
galaxy (single burst - filled red curve), SO galaxy (r=lGyr - 
orange dotted line), and late type galaxy (r=7Gyr - blue dotted 
line) are plotted for comparison. A young starburst (assuming 
L v oc - cyan dotted line) is also plotted. Different morpho- 
logical classes are represented by blue circles for spiral disks, 
green squares for peculiars, black diamonds for compacts, and 
red triangles for mergers. The mean error for B435-Z850 plot 
is 0.03mag (O.Olmag for V606-Isi4 respectively), much smaller 
than the symbol size. 



ent morphological classes or between different redshifts can be 
identified. This suggests that our objects are mainly composed 
of a mix of stellar populations including old, intermediate, and 
young populations. This is indeed confirmed by spectroscopic 
studies (see e.g. Hammer et al. 2005). This mix clearly appears 
on the color maps where the resolved colors allow us to inves- 
tigate the distribution of stellar populations within each object 
(see Fig.[T2b. As described above, the spiral disks share an old 
and red central population with a younger and blue population 
located in the disk. In the following, we try to characterize these 
populations quantitatively. 

6.3.1. Central color 

To derive the central color of the rotating spiral disks, we fol- 
low the procedure defined in Ellis et al. (2001). We construct a 
large elliptical aperture using second-order moments, and then 
a circular aperture with a radius of 5% of the semi-major axis 
of the large aperture is used to integrate the inner colors. Ellis 
et al. (2001) demonstrate that, following this method, the con- 
tamination from disk is negligible, and then, that these inner 
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colors can be associated with the color of the bulges. We show 
these 'bulge' -integrated colors as well as the modeled evolu- 
tion curves in Fig. [9] (see red filled dots). As an estimation of 
the uniformity of the color inside each bulge, we overplot the 
maximal pixel-by-pixel range of colors. 
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Fig. 9. B435-Z850 (top) and V606-Isi4 (bottom) observed color of 
the bulges for the rotating spiral disks sub-sample. The different 
evolutionary models are the same as in Fig. [8] Error bars are 
used to show the maximal pixel-by-pixel color range within 
each bulge. 

The modeled evolution curves are mainly used to com- 
pare the bulge integrated colors with the passively evolving 
system formed in a single burst at z-5 (elliptical model). 
Indeed, bulges as red as the elliptical model indicate that their 
populations have comparatively similar ages and metallicities 
(Pel etier et al. 19991 , and then that these bulges were formed 
at high redshift, whereas bluer bulges probably host a younger 
population. Moreover, dust extinction may significantly redden 
the color, and the galaxy sample certainly suffers from redden- 
ing effects because the mean inclination is particularly high ( 
~ 60°). This would make the bulges intrinsically bluer than 
observed. Nevertheless, even with the effect of dust extinction, 
we find that only 3 bulges are as red as the elliptical model, 
whereas most of them lie near the Sbc model evolution curve. 
This clearly suggests that these bulges host a significant popu- 
lation of young/intermediate age stars. 

6.3.2. Color of disks 

The disks of these distant rotating spirals show very blue out- 
skirts with the color consistent with young star-forming regions 
(see Fig. \12\i . Accounting for the uncertainties, we identify as 



a pure starburst region all those pixels where the observed B-z 
color (observed V-I ) is lower than lmag (0.4 mag) for CDFS 
(respectively HDFS and CFRS) galaxies (see Fig. [8]). By do- 
ing this, we are able to (i) compute the fraction of light coming 
from pure starburst regions, and (ii) find where this starburst 
activity is located in the disk. 

To derive the relative strength of the starburst activity, we 
first computed the ratio between the UV (observed B or V 
band) flux coming from the pure starburst regions over the to- 
tal UV flux. The total UV flux were defined to include all the 
pixels that are Act above the background, following the method 
described in Zheng et al. (2004). We find that between 5% and 
40% of the global UV flux can be associated with pure star- 
forming regions (mean value = 15%). 
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Fig. 10. Top: Three color images of three rotating spiral disks 
for which we have identified the star-bursting regions (B-z < 1 
or V-I<0.4). Middle: Color profiles for the CDFS galaxies. 
Bottom: Repartition of the starburst light for a growing ellipti- 
cal aperture. For each aperture, it represents the integrated UV 
light coming from pure starburst regions normalized by the UV 
flux coming from all star-bursting regions. 

In Fig. [10] (top) we show a three-color image of three 
galaxies for which we overplot a contour around the starburst- 
ing regions. It clearly appears that the starbursting regions 
are located towards the disk outskirts. More precisely, in Fig. 
[TOl (bottom), we define several elliptical apertures growing 
from the center of each object to an outer radius. The maximal 
radius being defined to include all the galaxy pixels that are Act 
above the background. This allows us to define a normalized 
radius and compare together galaxies with different sizes. 
Notice that this maximal radius roughly corresponds to the 
optical radius, as we find it to be ~ 3.2 times the disk scale 
length dPersic et al. 199 it . In each aperture, we integrate the 
UV flux coming from the pure starburst regions and normalize 
this value by the UV flux coming from all star-bursting regions 
(see above). These results are shown in Fig. QJJ (bottom) for 



12 



B. Neichel et al.: IMAGES II. Morpho-kinematic 



each galaxy (dotted lines), as well as the mean trend (solid 
line). This plot also clearly indicates that the starburst activity 
is increasing towards the disk outskirts. Assuming the maximal 
radius to be 3.2 times the disk radius, we find that more than 
95% (respectively 45%) of the starburst activity is located at a 
greater radial distance than Rd (respectively 2Rd). 



6.4. Stellar mass and SFR 

Table [3] shows the stellar mass and SFR for the rotating spi- 
ral disks sample. The stellar masses have been derived by 
Ravikumar et al. (2007) for galaxies located in the CDFS and 
by Flores et al. (2006) for the galaxies in the CFRS and HDFS. 
For the SFR, we report the SFRfjy and SFR/s from Hammer 
et al. (2005) for the CFRS galaxies. SFRjyy are derived from 
H„ luminosities and SFR/« from ISOCAM observations. For 
the CDFS galaxies, we derived the SFRrjy using the Kennicutt 
(1998) calibration based on the rest frame 2800 A and the 
SFR/« from Spitzer observations. To derive the SFR/« we used 
the 24yum MIPS data. At z ~0.6, this flux roughly corresponds 
to the rest-frame 15/im flux so we can use the Chary & Elbaz 
calibration to derive the total IR luminosity, and then use the 
classical Kennicutt calibration to derive the SFR/s. Although 
the SFR for CFRS and CDFS galaxies are derived with differ- 
ent methods, they give consistent results. Moreover, the ma- 
jority of the galaxies come from the CDFS that have the best 
quality data. For the galaxies that have not been detected by 
Spitzer or ISOCAM, only an upper limit of the SFR/« can be 
derived from the detection limits. 

Finally, in the last column of Table [3] we report the ratio of the 
stellar mass to the SFR ro/ , where SFR /o; =SFR/r+SFR[/i/ (and 
only a lower limit for galaxies not detected in IR). This ra- 
tio provides a time-scale that can be roughly associated with a 
mass-doubling time (if SFR is constant). 

In the previous section, we have shown that starburst re- 
gions are mostly located in the disk outskirts, suggesting very 
active and recent star formation in the outer parts of the disks. 
In a crude analysis, we can compare the M jr( ,#/SFR ratio be- 
tween these outer regions and the inner parts of the disks. We 
restrict this analysis for CDFS galaxies that have been detected 
by Spitzer, because they represent the most robust available 
data. Our purpose is not to quantify these values precisely, but 
rather to derive general trends. We first assume that the spa- 
tial distribution of the SFR2800A follows the observed B-band 
flux distribution: for z ~0.6 galaxies, the observed B-band in- 
deed roughly corresponds to the 2800A flux. To do so, some 
cruder assumptions are required, such as (i) the SFR/r follows 
the same distribution as the SFR2800A, and (ii) the stellar mass 
is sampled by the observed z-band flux. These last two as- 
sumptions are probably very approximate, since we have no 
idea how the dust is distributed within the galaxies. However, 
these assumptions are needed to distribute the SFR and the stel- 
lar mass at different radii and to evaluate the distribution of 
the M srf ,/;/SFR time scale. We computed this time-scale in dif- 
ferent annuli, from the center to the maximal radius. Results 
are shown in Fig. QT| All the galaxies show a negative gra- 



dient from large time scale near the center (between 4 and 8 
Gyr) to shorter scales near the disk outskirts (between 1 and 2 
Gyr). In the regions even farther out, the stellar populations are 
dominated by young starbursts (see Fig. ITOb . For these short- 
lived young stellar populations (typically <100 Myrs), a signif- 
icant correction of the mass estimation should be applied (see 
IBell et al. 2003b . Rather than try to correct for this effect, we 
arbitrarily decrease the M Jfc ;;/SFR time scale toward 0. 1 Gyr at 
the larger radius. 

Under all the assumptions made above, these observations 
clearly suggest an inside-out formation process in which the 
disk is gradually built from the outskirts. 
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Fig. 11. Ratio between stellar mass and SFR,„, as a function 
of galaxy radius. SFR, , is defined as the sum of SFRf/y and 
SFR/r. Only galaxies which have been detected by Spitzer are 
used in this analysis. The vertical line is at 1/3.2 of Rmax which 
roughly corresponds to the disk scale length limit. 



7. Discussion 

7.1. A low fraction of rotating spiral disks atz -0.6 

It is tempting to define morpho-kinematical classes from the 
relatively good coincidence between our morphological classi- 
fication and the kinematic classes (see Fig. [5]). Most galaxies 
can be distributed between three main distinct classes: (i) the 
rotating spiral disks for which both kinematics and morphol- 
ogy are relaxed, (ii) the galaxies with irregular morphology and 
perturbed kinematics, and (iii) the fully unrelaxed objects for 
which both morphology and kinematics are complex. 

It is useful to compare this distribution with what 
we know from the local universe. There, the overwhelm- 
ing majority of objects are dominated by rotation or 
dispersion (e.g. GASPH: IGarrido et al. 2002I - SAURON: 
Emsellem et al. 2004K whereas galaxies with anomalous kine- 
matics are generally ones in close pairs or interactions (e.g. 
Kannappan & Barton 2004). Disk galaxies constitute the ma- 
jority of local objects with ~ 70% of the intermediate-mass 
galaxies (Nakamura et al. 2004I IHammer et al. 2007I ) and we 
can reasonably assume that these disks are almost all dynami- 
cally relaxed as they lie on the tight Tully-Fisher relationship. 
Indeed, Kannappan et al. (2004) show that, at low z, large TF 
offsets are produced by kinematically anomalous objects. At 
higher redshift, Paper III find a similar result as they observe 
that the dispersion around the TF relation increases from RDs 
(0.31 mag) to PRs (0.80 mag) and CKs (2.08 mag). Restricting 
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ID 




SHR^y y 


SFR,« 


M„./;/SFR 










Gyr 


J033212.39-274353.6 


10.61 


1.0 


16.2 


2.4 


J033219.68-275023.6 


10.88 


5.1 


18.6 


3.2 


J033230.78-275455.0 


10.66 


2.9 


<10 


>3.5 


J033231.58-274121.6 


10.16 


3.2 


<10.5 


> 1 


J033237.54-274838.9 


10.70 


7.9 


22.3 


1.7 


J033238.60-274631.4 


10.53 


2.1 


<7.5 


>3.5 


J033226.23-274222.8 


10.72 


3.2 


17.5 


2.5 


J033248.28-275028.9 


10.09 


2.0 


<3.5 


>2.2 


HDFS4020 


9.82 








CFRS030046 


10.64 


3.9 


<15 


>2.3 


CFRS030085 


10.21 


2.3 


50.9 


0.3 


CFRS030619 


10.63 


1.6 


<13 


>2.9 


CFRS031353 


10.85 


4.1 


<25 


>2.4 


CFRS039003 




12.12 


94.95 





the distant TF relation to dynamically well-relaxed RDs, the 
local and distant relations have comparable dispersion (e.g. 



Pizagno et al. 2007 for the local reference). It suggests that 
both distant RD and local spirals are comparably kinematically 
relaxed. Of course, some kinematically distorted or even com- 
plex velocity fields may also fall on the TFR (see for instance 
Fig. 8/9 in Flores et al. 2006), but they represent a very small 
fraction. Then, a significant change between distant and local 
galaxies concerns this fraction of rotating spirals. Indeed, we 
found that at z ~0.6, ~ 1/3 of galaxies are regular rotating disks 
(see Section |2] and Tabled, which implies that the number of 
rotating spiral disks has increased by a factor ~2 during the last 
6 Gyrs. 

However, one could argue that our classification scheme 
(either morphological and kinematical) is too restrictive and 
that perturbed rotating Pec/Irr galaxies should be included 
in the rotating spiral disk sample. Even if we include this 
category among the rotating spiral disks, this would lead to a 
fraction of ~3/7 of "regular galaxies" that is still less than the 
local fraction. Moreover, Paper I convincingly showed that the 
kinematic classification is robust and that discrepancies seen in 
velocity fields or sigma maps correspond to the morphological 
perturbation seen in the optical images. Given the good 
correlation between kinematics and morphology, it is unlikely 
that some rotating spiral disks may have been missed. Only 
two galaxies that are morphologically classified as spiral disks 
are PR (J033226.23-274222.8 and J033248.28-275028.9) 
and one is CK (J033213.06-274204.8 - see also Sect. 4). 
These particular cases point out that morphology alone is 
insufficient for detecting a slight kinematic peculiar structure, 
such as is caused by minor mergers (e.g. Puech et al. 2007b 
for J033226.23-274222.8). On the other hand, over the 15 
RD galaxies, three are not morphologically classified as spiral 
disks (J033234.04-275009.7 ("Peculiar - Irregular" object in 
Fig. |U>, J033241. 88-274853.9 (first object in the third column 
of Fig.0, J033245.1 1-274724.0 a compact galaxy). 
These exceptions show that some perturbed morphology 
can nonetheless have relaxed dynamics and may require 



further analysis (e.g. Puech et al. 2008b for a detailed anal- 
ysis of J033241. 88-274853.9). In each instance, to verify 
whether k-correction may have influenced our classification, 
we re-examined the morphology on the sole basis of the 
F850LP/F814W filter. Doing so, we find that all the morpho- 
logical disturbed features that prompt us to classify an object 
as peculiar (and which are generally blue) are always present 
in these 'red' filters. This suggests that k-correction does not 
strongly affect our conclusions. 

Finally, the main limitation of this study is the still small 
number of objects, leading to possible non negligible statistical 
uncertainties. However, each of the 4 different fields taken 
independently give the same distribution in morphological as 
in kinematical numbers. A simple^ 2 test between the different 
fields confirms this trend. Moreover, the sample follows the 
luminosity function in the z = 0.4 - 0.75 range (see Fig. 1 in 
Paper I). This indicates that, even if the number of objects is 
small, our result, a small number of rotating spiral disks at z= 
0.6, is robust. 



7.2. Comparison with previous works 

The evolution in the number fraction of regular spirals is sur- 



prising and in contradicts former studies (e.g. Lilly et al. 1998 



: IRavindranath et al. 20041 |Sargent et al. 2007| > that found that 
most of the spirals were in place at z=l. The differences ob- 
served with these former studies certainly come from the dif- 
ferent methodology used to identify disk galaxies. 

Lilly et al. (1998) used a 2D light-fitting procedure, approx- 
imatively similar to the one presented in Sect. 13.1.21 The main 
differences are that (i) they used a de Vaucouleur light profile 
for the bulge components instead of a Sersic profile, and (ii) for 
each galaxy, they subtracted a version of itself rotated by 180° 
from it. This last step is performed to "symmetrize" the galax- 
ies, but it would not modify the global underlying light distri- 
bution. They then classified as disk-dominated objects all those 
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galaxies with B/T<0.5. Even if their method is slightly differ- 
ent from our derivation of B/T, we have seen in Sect. [5] and in 
Fig.[6]that a method only based on this structural parameter is 
not able to clearly distinguish RD from CK so it unavoidably 
overestimates the number of spirals. 

In their analysis, Ravindranath et al. (2004) use a single 
Sersic component to model the brightness profile of the galax- 
ies. They then classify as disk-dominated all objects with n<2. 
In a crude approach, we reproduced their procedure for the 
galaxies located in the CDFS, following the method described 
in Sect. 13.1.21 but with only one Sersic component per galaxy. 
By doing this, we find 30 galaxies with n<2, over which 8 are 
RD but 13 are CK (9 are PR). This is agree with Cassata et 
al. (2005), who find that "the sersic index alone provides just 
a broad and not unequivocal indication of the morphological 
type". 

The method developed in the present study is not fully auto- 
matic but instead based on a detailed analysis of each galaxy in 
order to assign them a morphological label. Structural param- 
eters (Ri w if and B/T), color information, and visual inspection 
are used in a complementary way. In this process, only the more 
regular galaxies (in their morphology and color distribution) 
are identified as spirals. This method indeed leads to a smaller 
number of spirals compared to previous studies; however, the 
classification seems robust because, on one hand, almost all 
these spirals are rotating disks, and on the other, almost all the 
rotating disks are spirals. This indicates that our classification 
is, at least, efficient in correctly recovering the rotating spirals. 
We then suspect that, among the objects identified as spirals by 
former studies, a large fraction of them are dynamically per- 
turbed or complex and have indeed peculiar morphologies. As 
an illustration of this, 12% of the objects classified as large 
disks galaxies by Lilly et al. (1998) are included in the list of 
visually identified mergers presented by Le Fevre et al. (2000). 
As discussed in Sect. [5] automatic methods unavoidably mix 
rotating spirals with dynamically complex objects. Strikingly, 
we precisely find that only about 50% of the objects identified 
as spirals by these methods are in fact rotating spiral disks. It is 
then probable that former studies have overestimated the num- 
ber of spirals, even if only a braoder 3D spectroscopic survey 
could definitely resolve this issue. 

7.3. Are rotating spiral disks forming stars inside out? 

In this section we consider the subsample of rotating spiral 
disks with emission lines, e.g. 16% of the whole sample of 
M sle u > 1.5xlO 1() M , galaxies at z=0.4-0.8. Because it is prob- 
able that these galaxies are evolving into present-day spirals, 
they vary likely represent one fourth of today's spiral progen- 
itors (see Sect. 17.11 ). Based on their morphology, kinematics, 
and color distribution, the rotating spiral disks identified in this 
paper share many properties with local spirals. They all show 
regular patterns, with a red central bulge surrounded by a blue 
disk, and all have integrated colors between the late type and 
Sbc models. 

However, those disks are certainly not passively evolving 
systems because they are actively forming stars (Table [3]), and 



their bulges are bluer by ~1 magnitude than elliptical stellar 
population models (Fig. [9]). The latter result is not new, since 
its discovery by Ellis et al. (2001), who concluded that sec- 
ondary star formation superimposed on preexisting old popula- 
tions (rejuvenation) have occurred in the inner regions of spi- 
rals galaxies (Thomas et al. 2006). Alternatively, Hammer et 
al. (2005) suggest that bulges at z ~0.6 may have been formed 
by mergers, few Gyrs prior to the redshift epoch. 

Let us assume that the subsample of rotating spiral disks 
is made of galaxies that are likely unaffected by galaxy inter- 
actions, at least since z=l. It is quite likely that they have not 
undergone a major merger recently, and the probability that a 
galaxy at z=0.6 experiences another merger event is quite low. 
Table [3] shows that most of them (if not all) present a doubling 
time ranging from 0.2 to 4 Gyrs. In the absence of interactions, 
it may be reasonable to assume an exponential decay of their 
SFRs, or a constant SFR as a lower limit of their past SFR. 
Because all of them but one are disk-dominated, it means that 
most of the stars in their disks have been formed rapidly, e.g. in 
a few Gyrs or less. Thus it implies quite a rapid building of the 
stellar disk and, also assuming a pre-existing old population in 
the bulge, a significant decrease in their B/T values. 

What could be the subsequent evolution of the rotating spi- 
ral disks at z=0.4-0.8? Sect. l6.4l provides evidence that recent 
starbursts are occurring in their outskirts. We also see decreas- 
ing mass-doubling time from the center to the edge, which sug- 
gests an inside-out disk growth; however, this conclusion is 
severely dependent on the assumption that IR light follows the 
UV light, which is difficult to verify given the limited spatial 
resolution of Spitzer. An independent argument is provided by 
their location along the Tully Fisher relation (TFR, see Paper 
III). They delineate a TFR that is 0.34 dex shifted towards 
lower stellar masses than are local spirals. This is in remarkable 
agreement with the evolution of the stellar mass metallicity (0.3 



dex on M-Z, Liang et al. 2006[ ). It may mean that the rotating 
spiral disks have to double their stellar masses to reach the local 
TFR and the M-Z relationship. Combined with the rapid forma- 
tion of their stellar disks prior to their redshift epochs and with 
the fact that they show recent starbursts in their outskirts, it 
seems likely that their stellar disks are still growing inside-out 
(and their B/T values are decaying) during the 6 Gyr elapsed 
time since z=0.6. Such a conclusion is also in good agreement 
with results from Trujillo et al. (2005). 

7.4. Testing various scenarios of spiral formation 

As these galaxies are dynamically relaxed, we expect the SFR 
to be self-regulated by the disk (e.g. Silk 2003) and star forma- 
tion to decay monotonically as the gas supply is exhausted. In 
contrast with theoretical predictions, we find that one quarter 
of the local spiral progenitors are sustaining a high SFR, which 
necessarily implies massive gas accretion or, alternatively, pre- 
existing large amounts of gas in an already relaxed gaseous 
disk. What is the origin of this gas and which scenario of spiral 
formation is able to reproduce our observations? 

Despite the still small number of observations, these distant 
disks seem to have systematically lower V/<x values than their 
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local counterparts (Puec h et al. 2007l >. They seem to be heated, 
which could possibly be associated with turbulence induced by 
processes related to star formation, and unlikely to minor merg- 
ers given that our observations are quite good at identifying mi- 
nor mergers (Puec rTet al. 2007bt . At the same time that rotating 
disk spirals are transforming their gas to stars to reach the local 
relationships (TFR and M-Z), they have to evolve towards thin 
disks, e.g. higher values of V/cr. 

We may infer that the necessary gas that feeds the star for- 
mation has a specific orientation to stabilize their disks dur- 
ing the past 6 Gyrs. Processes based on isotropic gas accre- 
tion (Dekel & Birnboim 2006 ) may have difficulty reproducing 
such an evolution, while gas accretion from clumps or filaments 
(e.g. Keres et al. 20051 1 requires some preferential alignment to 
stabilize the disk. Thus it may be necessary for the gas to al- 
ready pre-exist in a gaseous thin disk, and we may be simply 
seeing the gas transformation into stars in such systems. 

An alternative is that this gas reservoir is provided by 
re-accretion following a past major merger (Barnes 2002; 
Hammer et al. 2005IGovernato et al.~20 07 ). This post-merger 
gas accretion will constantly deposit gas at the outer edge of a 
spheroidal remnant, which will naturally give rise to the inside- 
out formation of a new disk surrounding the newly formed 
bulge. Extrapolating Barnes (2002) simulations at higher red- 
shift, where the gas fraction is expected to be higher, may re- 
sult in remnants with even larger and more massive gas disks. 
Following this scenario, the morphological expectations are a 
relatively undisturbed disk (since it is assembled progressively 
by smooth re-accretion), with a clear inside-out signature with 
the youngest stars located at the outer regions of the disk, which 
closely matches what is observed in the rotating spiral disks 
(see above). From the viewpoint of the kinematics, we expect a 
relatively high velocity dispersion of the gas caused by shock- 
heating between the infalling gas material with the already 
accreted gaseous disk, and a progressive smooth increase in 
the angular momentum of the disk with time. Strikingly, these 
kinematical properties have already been observed in these 
galaxies by Puech et al. (2007). Besides this, the bulge colors 
shown in Fig. |9]are also easy to reproduce if these bulges have 
been assembled a few Gyrs prior to the epoch at which we are 
observing them. 

Since the disk growth time scales range from less than 1 
Gyr to 4 Gyr, this suggests that these rotating spiral disks would 
have undergone their last major merger at z=l or slightly later. 
At these redshifts, the gas fraction in galaxies is expected to 
be much higher (Liang et al . 2006| ILaw et al. 20 07 ), which is 
more than two times higher at z=0.65 than today. This gives 
even more credit to the post-merger gas accretion scenario, be- 
cause the gas reservoir available for the disk rebuilding phase 
would be fairly large. From z-1 to z-Q, the merger rate in 
this mass range is expected to be between 50 to 75% (see 
introduction in IHammer et al. 2007i l. Then, between 25 and 
50% of local intermediate-mass galaxies would have escaped 
such an event since z=l, which is, given the large uncertain- 
ties associated, in relatively good agreement with the estima- 
tion of the fraction of intermediate-mass rotating spiral disks, 
i.e., 33%. Given the much higher merger rate expected at z >1 
dConselice et al. 20031 ILotz et al. 2006i l, our results are reason- 



ably consistent with these spiral disks possibly undergoing a 
major merger before z- 1 . 

8. Conclusion 

We have presented a morphological analysis of a representative 
sample of intermediate-mass galaxies at z ~0.6. To derive our 
morphological classification, both structural parameters and 
color maps were used, combined with a visual inspection of 
each object. The method used in the present study is slightly 
different from previous studies, as we made used of the 
color information to distinguish only those spiral disks with 
properties comparable to local spirals (bulge redder than disk). 
This method was tested against an independent kinematical 
classification, and we find good agreement between the 
morphological and the dynamical states of the galaxies (see 
Fig. |5J. There is a definite offset in morphology between RD 
and CK. The major overlap between the different kinematical 
classes are for irregular galaxies known to span a wide range 
of kinematical properties. 

We also computed the concentration, asymmetry, GINI and 
M20 parameters and find they are not able to recover the same 
correlation. It appears that, when applied automatically, these 
methods will miss around half of the complex kinematics 
and unavoidably overestimate the fraction of spirals. These 
results show that, when derived properly, the morphological 
information can be representative of the underlying kinemat- 
ical properties, even at z ~0.6. This correlation represents an 
interesting new tool for understanding the mechanisms in the 
formation and evolution of galaxies. 

Our classification scheme allows us to isolate a particular 
morpho-kinematical class: the rotating spiral disks. This class 
includes objects sharing the spiral disk morphology and rotat- 
ing disks properties. We find that the number of these rotat- 
ing spiral disks is much lower (by a factor of ~2) at z ~0-6 
than today. The radial distribution of color along their disks 
and their high SFR suggest that their stellar disks have been 
formed relatively recently before z ~0.6 through inside-out 
processes. These findings are consistent with scenarios for 
which the large gas supply stabilizes the disk, e.g. leading ul- 
timately to present-day thin disks by increasing V/cr values 
with time. It opens up two possibilities, either that z=0.6 ro- 
tating spiral disks are forming their stars from a pre-existing 
gas reservoir already distributed in a thin disk or, alternatively, 
a post-merger accretion scenario where a disk is rebuilt thanks 
to gas accretion left over from the merging event (Barnes 2002 
; IHammer et al. 2005l l. 
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Table 4. Catalog of derived parameters for the working sample. 
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a Redshift measured by [Oil] emission. 
* Bulge fraction derived from Galfit. 

c Quality factor: 1 = secure; 2 = possibly secure; 3 = fit failed or unreliable 

d Morphological galaxy type: Sp = Spiral Disk; Pec/Irr = Peculiar-Iregular; Pec/T = Peculiar-Tadpole 

Pec/M = Peculiar-Merger; C = Compact; M = Merger. 

' half-light Radius (kpc). 

J Disk radius (kpc) derived from Galfit. 

g Kinelatical type: RD = Rotating Disk; PR = Perturbed rotation; CK = Complex kinematics. 
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Fig. 12. F850LP/F814W image, color map. Explanations are given at the end of this figure. 
Rotating spiral disks are marked with an asterisk near their name. 
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Fig. 12. F850LP/F814W image, color map. Explanations are given at the end of this figure. 
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Fig. 12. F850LP/F814W image, color map. The size of each image is 20x20 kpc. Galaxies observed with ACS have F850LP 
images, galaxies in CFRS and HDFS have F814W images. For each target, the name and redshift are labeled at the top-left 
and top right in the F850LP/F814W image. Rotating spiral disks are marked with an asterisk near their names. Morphological 
and kinematic classifications are indicated at the bottom-left and bottom right, respectively. Color maps are B435-Z850 for CDFS 
objects and V 60 6-l8i4 for objects in HDFS and CFRS (except CFRS039003, CFRS030488 with B 4 g ( )-I 81 4 and CFRS031032 with 
V535-I814). All the galaxies are aligned with the GIRAFFE IFU (see Paper I). 



